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Abstract 

We present all relevant details of our calculation of the complete next-to- leading order 
{0(a 2 s a)) QCD corrections to heavy flavor photoproduction with longitudinally polarized 
point-like photons and hadrons. In particular we provide analytical results for the virtual 
plus soft gluon cross section. We carefully address the relevance of remaining theoretical 
uncertainties by varying, for instance, the factorization and renormalization scales inde- 
pendently. Such studies are of importance for a meaningful first direct determination of 
the polarized gluon density Ag from the total charm production spin asymmetry by the 
upcoming COMPASS experiment. It is shown that the scale uncertainty is considerably 
reduced in next-to-leading order, but the dependence on the charm quark mass is sizable 
at fixed target energies. Finally, we study several differential single-inclusive heavy quark 
distributions and, for the polarized HERA option, the total bottom spin asymmetry. 



1 Introduction 



Measuring the unpolarized gluon density of the nucleon g(x, /i 2 ) at a scale \i as a function 
of the momentum fraction x presents considerable theoretical and experimental challenges 
and thus serves as a benchmark for the steady progress in QCD. The determination of 
g(x, /i 2 ) from measurements of the structure function F 2 in deep inelastic scattering (DIS) 
is hampered by the absence of direct couplings to the electroweak probes (7*, Z, W^). 
However, the increasingly precise F2 data from HERA |I| still serve to constrain the small- 
x behaviour of g(x, // 2 ) indirectly in the region 10~ 4 < x < 10~ 2 with an accuracy of about 
10% p| from the observed scaling violations dF 2 (x, /i 2 )/<9/i 2 . To determine g(x,fi 2 ) over 
the entire x region, i.e., also at larger values of x, studies of exclusive reactions like direct 
photon or di-jet production, where the gluon already enters in leading order (LO), are 
indispensable. Such measurements are often experimentally much more involved and less 
precise than inclusive DIS. Nevertheless, our knowledge of the unpolarized gluon density 
has greatly improved in the past few years (for recent QCD analyses see [0]), but in the 
region x > 0.1 the situation is still far from being satisfactory. Here the uncertainty in 
g(x,fi 2 ) easily amounts to about 100% @, ||, [J- 

Concerning the spin properties of the gluons in a longitudinally polarized nucleon, 
the unpolarized gluon density g(x,n 2 ) is defined as the sum of the two possible helicity 
distributions, whereas the corresponding polarized gluon density Ag(x,fi 2 ) is given by the 
difference. In general we have for a parton / with / = g,q,q 

unpolarized : f(x,/j, 2 ) = / + (x,/i 2 ) + /_(a;,/i 2 ) , 

polarized: A/(x,/i 2 ) = f + (x, /i 2 ) - fi 2 ) . (1) 

Here /+ and /_ are the densities with the parton spin aligned and anti-aligned to the spin 
of the nucleon, respectively. In order to measure the two independent combinations in ([!]), 
we need experimental data and theoretical calculations distinguishing between different 
initial helicity states. 

The long list of spin-dependent DIS experiments f| and the recently completed next- 
to-leading order (NLO) framework for the evolution of the Af || I7|] may lead to the 
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expectation that the polarized gluon distribution Ag(x,fi 2 ) should be known with almost 
similar accuracy as g(x,fi 2 ) by now. This is, however, not the case as was revealed by 
all NLO analyses || || HH, °f presently available spin-dependent DIS data. In fact 
it turned out that the x-shape of Ag is even almost completely unconstrained. This 
ignorance is, of course, also reflected in present values for the first moment of Ag(x,fi 2 ), 
defined by 



which can be estimated at best with an error of 100% for the time being. Ag(fi 2 ) plays 
an important role in our understanding of the spin- 1/2 sum rule for nucleons 



where AE is the total polarization carried by the quarks and antiquarks and L z denotes 
the sum of the non-perturbative angular momenta of all partons. 

There are three main reasons for the present problems to pin down Ag(x,fi 2 ): 

• The measurements of the nucleon spin structure function gi, the polarized analogue 
to the unpolarized structure function Fx, are still in a "pre-HERA" phase. The 
kinematical coverage of the fixed target experiments [f| is by far not sufficient to 
constrain Ag(x,[i 2 ) from scaling violations dgi(x, ^L 2 )/dfi 2 . 

• As already mentioned, the unpolarized gluon density is also constrained by several 
exclusive reactions, but corresponding measurements in the polarized case are still 
missing. 

• A momentum sum rule for spin-dependent parton densities is lacking, i.e., we can- 
not infer any constraint on Ag from the already somewhat more precisely known 
polarized quark distributions. In addition, the spin- dependent parton densities Af 
defined in (P are not required to be positive definite. 




(2) 




(3) 



Nothing can be done about the last point, of course. The small- x region of g\ could be 
explored at HERA in case that the option to longitudinally polarize also the proton beam 
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fT2|j will be realized in the future. First measurements of Ag in exclusive reactions will 



be provided by the COMPASS fixed target experiment at CERN [13] and the BNL RHIC 



polarized pp collider [14| , which are both currently under construction. 



For the determination of the gluon distribution, heavy quark (Q = c, b) photoproduc- 
tion 

19 ^QQ (4) 

is an obvious choice (an arrow denotes a longitudinally polarized particle from now on). 
The reconstruction of an open heavy quark state is experimentally feasible, and in LO 
only the photon-gluon fusion (PGF) process in ([|) contributes, which may lead to the 
hope that an unambiguous determination of Ag can be performed. Thus open charm 
photoproduction will be used by the upcoming COMPASS experiment [O to measure 



Ag. All theoretical studies of (|j) have been performed only in LO so far [|T^, T7, 
However, LO estimates usually suffer from a strong dependence on the a priori unknown 
factorization and renormalization scales. Also there are new NLO subprocesses induced 
by a light quark replacing the gluon in the initial stateQ. Here the question arises if the 
PGF contribution (|]) still dominates in the experimentally relevant kinematical region 
as is desirable for a precise determination of Ag. Finally, the NLO corrections have been 
shown to be sizable near threshold in the unpolarized case Jl8|, |19[| . Clearly, a NLO 
calculation also for the spin-dependent case is warranted in order to provide a meaningful 
interpretation of the forthcoming experimental results. 

This paper provides all relevant details of the first calculation of the complete NLO 
(0(a 2 s a)) QCD corrections to heavy flavor photoproduction with point-like photons |20| . 



In |2(J we only highlighted some of the most important phenomenological aspects, but we 
skipped most calculational details. In addition, we now present, again for the first time, 
NLO studies of differential single-inclusive heavy quark distributions. In the next section 
we will first make some general technical remarks concerning the polarized calculation. In 
Section 3 we recall the known LO results and extend them to n dimensions as is required in 



furthermore, the on-shell photons in (|J) cannot only interact directly, but also via their partonic 
structure. However, LO estimates of this unknown "resolved" contribution are small for all experimentally 
relevant purposes mM. 
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course of the NLO calculation. In Section 4 we calculate the virtual one-loop corrections 
to dJ) and examine the gluon bremsstrahlung process —>■ QQg in detail. Section 5 
is devoted to the new genuine NLO contribution with a light quark in the initial state 
7<f — > QQq. The relevant formulae for calculating total and differential single-inclusive 
heavy quark cross sections can be found in Section 6, where we also present some further 
phenomenological studies. Finally, our main results are summarized in Section 7. In 
Appendix A we present the details of the phase space calculation. In particular, we focus 
on peculiarities which arise in a polarized calculation using dimensional regularization. 
Here we also supply the parametrizations of the parton momenta used in our calculations. 
Appendix B contains several helpful remarks concerning the calculation of the tensor 
integrals needed for the virtual corrections and Appendix C collects the analytical results 
for the polarized virtual plus soft cross section. 



2 Some General Technical Remarks 

In the calculation of the NLO corrections we will encounter the usual array of ultraviolet 
(UV), infrared (IR) and mass/collinear (M) singularities. We choose the framework of n- 
dimensional regularization to deal with all of these various types of singularities. Since our 



calculations proceed along similar lines as in |T8| for the corresponding unpolarized case, 
we adopt their notation n = 4 + e for the deviation from four space-time dimensions in 
order to facilitate comparisons of the intermediate results. Of course, all results presented 
here can be easily converted to the more common choice n = 4 — 2e by just replacing e — > 
— 2e accordingly. In the calculation we simply identify the dimension used to regularize 
the UV divergencies n < 4 with the one used for the IR divergencies n > 4. Thus we do 
not distinguish between Euv < and sir > 0, which for example leads to the following 
result for the basic loop integral: 

= for a > . (5) 

Choosing the n-dimensional regularization method introduces some complications when 
polarized processes are investigated due to the unavoidable presence of 75 and the totally 
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anti-symmetric Levi-Civita tensor e^ vpa . First we shall recall how these quantities appear 
when projecting onto the helicity states of the incoming particles, and then we will explain 
how to deal with them in n dimensions. 



One can calculate the squared matrix elements for both unpolarized and polarized pro- 
ses simultaneously using the squai 
h\ and h 2 of the incoming particles: 



cesses simultaneously using the squared matrix elements \M\ 2 (hi, h 2 ) for definite helicities 



unpolarized : \M\ 2 = - [\M\ 2 (++) + |M| 2 (+-)] , (6) 
polarized: A |M| 2 = l - [\M\ 2 (++) - |M| 2 (+-)] . (7) 

This is of course highly desirable, since we obtain an important consistency check by 



comparing with the already known unpolarized results [0, |I9 |. To obtain \M\ (hi, h 2 ) we 



use the standard helicity projection operators (see, e.g., |21| ) 

1 



e M (/ci,Ai)e*(A;i ) Ai) } 



h p h a 



-9tw + *Aie 



(8) 



for incoming photons with momentum k\ and helicity Ai (accordingly for gluons with k 2 
and A 2 ) and 

u(k 2 ,h)u(k 2 ,h) = ^ 2 (l-h l5 ) (9) 
for incoming quarks with momentum k 2 and helicity h (analogously for antiquarks). 

We note that in the unpolarized case one has to average over the n — 2 spin degrees of 
freedom for each incoming boson in n dimensions. This can be achieved by the replacement 
—1/2(^2, — > —l/(n — 2)g IJjV = —1/(2 + e) in (g) leaving (^j) unchanged. However, it 
is convenient, both for the calculation and for the presentation of the results, to define 
instead 

E ={ + ^ un P° larized / 1Q ) 
£ ~ 1 1 polarized 

One can then perform the unpolarized and polarized calculations using (Bp, if one multi- 
plies the results by a factor E e for each incoming boson. We have also always identified 
the additional four-vector if usually appearing in (pj) |21| to be that of the other incoming 
particle. This is possible since k\ ■ k 2 — s/2 ^ in (^) and simplifies the rather lengthy 
intermediate results considerably. 



p i-> o — — q i — ^ ^ _i_ ^ , ^ 



Figure 1: Graphical "rule" illustrating the replacement of the physical polarization sum 
(o) by —g^v (•) and appropriate ghost contributions (dashed lines). The minus signs in the 
lower half are due to the cut ghost loop. 

As a further simplification one can drop all terms other than — in the symmetric 
(unpolarized) part of (|8]). This of course means that unphysical polarizations will be 
kept in the polarization sums. However, unphysical photons decouple completely and 
unphysical gluons do not contribute as well, except for those subprocesses where one 
encounters a triple-gluon vertex. There one has to introduce incoming external ghost fields 
to cancel these unphysical parts p2| , when using —g^ instead of the physical polarization 
sum P pu = J2\ e fj.{^) e tW EH- Fig. [| illustrates this elimination of such terms by adding 
appropriate external ghost contributions. The extra factors (—1) multiplying each ghost 
contribution are due to the cut ghost loop. 

The quantities e pupa and 75 introduced by (^) and (^|), respectively, are of purely four- 
dimensional nature and there exists no straightforward continuation to n 7^ 4 dimensions. 



We treat them by applying the HVBM prescription [^3| which provides an internally 
consistent extension of e pvpa and 75 to arbitrary dimensions. In this scheme the e-tensor 
continues to be a genuinely four-dimensional object and 75 is defined as in four dimensions, 
implying {7^,75} = for /x = 0, 1, 2, 3 and [7^,75] = otherwise. This effectively splits 
the n-dimensional space into two subspaces, each one equipped with its own metric: 
one containing the four space-time dimensions and one containing the remaining n — 4 
dimensions, denoted "hat-space" henceforth. In the matrix elements we then encounter 
not only conventional n-dimensional scalar products of two momenta, like k-p = g^k^p", 
which can be expressed in terms of the usual Mandelstam variables, but also similar scalar 
products in the hat-space k ■ p = g^vk^p 11 ■ 
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These additional terms would complicate the phase space calculations considerably, 
but one can eliminate most of them by choosing the coordinate system appropriately. 
The three-body phase space integration for the NLO 2^3 processes exemplifies the 
problem: Since we are only interested in single-inclusive heavy quark cross sections in 
our calculation, it is possible to assign only the momenta of the two outgoing unobserved 
partons with non-vanishing hat-space contributions. The fully measured momenta of the 
initial states and of the observed heavy (anti)quark remain purely four-dimensional. A 
convenient choice of coordinates is the "Gottfried- Jackson frame" p4| , p5| , where the phase 
space integration is performed in the rest frame of these two unobserved particles. The 
remaining three momenta can then be oriented in such a way that they lie in, e.g., the y — z 
plane. To further facilitate the phase space integrations one can finally choose one of these 
three vectors to have only a non-vanishing ^-component (this freedom can be exploited 
to simplify the angular integrations). Using such a choice of coordinates, one thus ends 
up with only one scalar product of hat momenta k 2 , which simplifies the calculations 
considerably. In App. A one can find all required formulae concerning the phase space 
integration in the polarized case and the parametrizations of the parton momenta used in 
our calculations. For the additional integrals which depend on k 2 we furthermore show in 
App. A that thanks to the heavy quark mass they are all at least of 0(e) and hence drop 
out when the limit n — > 4 is taken in the end. Thus we arrive at the welcome conclusion 
that in the particular case of our calculation hat momenta terms in the matrix elements do 
not contribute to the final result and can be ignored^. Concerning all 2 — > 2 subprocesses 
presented in Sections 3 to 5 it is then obvious that the same holds true. Since three of 
four external particles have their momentum fully measured, the fourth is determined by 
energy-momentum conservation and thus all hat components can be eliminated from the 
calculation. 



2 This is different in calculations involving only massless particles, as was discussed in J26| in the 
context of direct photon production. Notice also that the unphysical helicity violation at the qqg- vertex 
in the HVBM scheme in n dimensions (see Q for details), is not relevant for our calculation, since this 
vertex does not appear in the mass factorization procedure (see Sections 4 and 5). 
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3 Born Cross Section in n Dimensions 



(a) 



(b) 



Figure 2: Feynman diagrams for the LO photon-gluon fusion process — > QQ. 



In this section we will recall the well-known LO results for the unpolarized and po- 
larized photoproduction of heavy flavors. Since we encounter 1/e 2 poles in our NLO 
calculation we have to extend these calculations up to 0(e 2 ) in n = 4 + e dimensions. For 
the contributions to the Born amplitude depicted in Fig. [|we use the following momentum 
assignment 

j(k 1 )+g{h)^Q(p 1 ) + Q(p 2 ) (11) 
and the corresponding Mandelstam variables are given by 



s = (h + k 2 ) 2 , 



U = t 



m 



(k 2 - p 2 ) 2 - m 2 



Hi 



u 



m 2 = (ki -p 2 ) 



m 



(12) 



where s + t x + u x = 0, k\ = ("on-shell" photon), k\ = 0, and p\ 



m 2 with m 



denoting the heavy quark mass. All trace calculations in n = 4 + e dimensions are per- 
formed using the package TRACER f2~7||. In order to present the unpolarized and polarized 



results simultaneously in the most compact form, we will use \M\ 2 here, and in the rest 



of the paper, to denote both the unpolarized \M\ and polarized A|M| color-averaged 
squared matrix elements calculated according to Eqs. (|B|)-(|lTi|). Similarly, in ( |T3| ) below 
Bqed denotes either the unpolarized Bq ED or the polarized ABq ED . The LO result can 
then be expressed as 



E 2 g 2 s e 2 e 2 Q B QED 



B, 



QED 



t\ Mi Am 2 s 

Mi t\ tiUi 



1 - 



m 2 s 



t t u x 



+ e 



Atxui 



(13) 



AB 



QED 



2m 2 s 

t X Ui 



where g s and e are the strong and electromagnetic coupling constants, respectively, and 
eg is the electromagnetic charge of the heavy quark in units of e, e.g., eg = e c = 2/3 for 
charm quarks. Notice that the polarized ABqed retains its four- dimensional form and 
receives no 0{e) contributions in contrast to the unpolarized Bq ED . 

Using the standard 2 — 2 phase space in n dimensions, 

dPS 2 = ^ [(4vr) 2+£ / 2 r(l + e/2)] ^ ^ ~ — ) ^ dfidui (14) 
one can then write the n-dimensional Born cross section as 

iJla ' ! ~ - F £ 5(s + t 1 + Ul )\M\ 2 , (15) 



dtidu\ 



2„\ £ / 2 



F E , J[(4.)^r(l + £ /2)]-(^^i) 

where _F £ collects all phase space factors given in (|I4]), the flux factor l/2s, and the 
mass parameter /x introduced to keep the gauge couplings g s and e dimensionless in n 
dimensions, a denotes the unpolarized and polarized cross section a and Act, respectively. 
If one is only interested in the Born result itself, one can of course perform the e — » 
limit in ( ]T3|) and simply use F £= o = l/(16irs 2 ). Our four- dimensional results for a$ and 
Acjgj (O) agree with those in [L8| and p^, p0[ , respectively. 



4 NLO Gluon Contribution 

Next we turn to the NLO corrections to the PGF process (Q), where one-loop virtual and 
gluon bremsstrahlung contributions have to be taken into account. The one- loop virtual 
corrections displayed in Fig. ^ have the same 2^2 kinematics as the Born graphs in 
Sec. 3 and can be also calculated using ([TT|) , (|12|) , (|14D , and fllSP . At 0(aa 2 s ) only the 
interference between the virtual (V) and Born (B) amplitudes of Figs. [2] and ^ contributes 

\M\ 2 VB = 2Re (MyM|) = E 2 g 4 s e 2 e 2 Q \lC F V QED + C A V OK ] , (16) 

where all quantities with a tilde denote again, as in Eq. (|13D , both the unpolarized and 
polarized expressions, e.g., Vqed denotes either Vqed or the spin-dependent AVqed- The 





Figure 3: The NLO virtual corrections to jg — > QQ. Reversing the heavy quark lines, except 
for the non-planar graph (b), yields the remaining graphs. Massless particle loops similar to 
graph (i) vanish, see App. B. 

color factors associated with the abelian and non-abelian parts are Cp = (Nq — 1)/ (2Nc) 
and Ca = Nq, respectively, where the number of colors is Nc = 3 for QCD. We note 
that Vqed, which receives contributions only from the graphs (a), (c), (d) and (f)-(h) in 
Fig. H], corresponds to the process where the gluon is replaced by a photon in the initial 
state, i.e., 77 — > QQ. A complete NLO QCD 0(a 2 ot s ) calculation of this process has 
been performed recently in |30|] for both the unpolarized and polarized case. Our NLO 
results for the QED-part of 7$ — > QQ agree analytically with the ones presented in [|30|j . 

In the loop-calculations we encounter Feynman integrals with up to four propagators 
in the denominator. We define the corresponding one-loop scalar one- (Aq), two- (Bo), 
three- (C ) and four-point (D ) functions as in Ref. e.g., the four-point function 
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needed for the box graphs in Figs. |3| (a), (b) is defined by 

D (qi, q 2 , <?3, m x , m 2 , m 3 , m 4 ) = (17) 
_ £ f d n q 1 
^ J (2ir) n (q 2 - m\) [(q + q{) 2 - mj] [(q + q 1 + q 2 ) 2 - mj] [(q + q x + q 2 + q 3 ) 2 - m 2 } ' 

where the four external momenta satisfy qi + q 2 + q 3 + q^ = and the m, are the internal 
masses. The required scalar integrals are conveniently collected in [EH], however, we have 
checked them using the standard Feynman parametrization techniques. Each fermion 
propagator and each triple-gluon vertex in the loop introduces a loop momentum q^ in the 
numerator. A glance at Fig. then reveals that the maximal number of loop momenta we 
face in the numerator is one less than the number of propagators, except for graph (i). In 
particular, one has to deal with tensorial four-point integrals of first to third (q^q u q p ) 
order and with tensorial three- and two-point integrals of first and second (q fl q u ) 
order. We have developed a program which automatically reduces these tensor integrals 
to a set of scalar ones by using an adapted Passarino-Veltman decomposition method 
|32|| , which properly accounts for all possible n-dimensionally regulated divergencies in 



QCD. Since this procedure is quite common, we will just mention a few helpful details in 
App. B. 

In the virtual cross section UV, IR and M singularities show up as 1/e poles. In the 
non-abelian Oi-T-part also double poles 1/e 2 occur when IR and M singularities coincide. 
The UV divergencies are removed by the renormalization procedure, which we implement 
using the common counterterm method ( "renormalized perturbation theory"). The coun- 
terterms introduce additional contributions similar to those in Figs. |^ (c)-(i), but with the 
loops replaced by corresponding renormalization constant dependent "interactions" . For 
the internal gluon propagators we use the Feynman gauge. As the renormalization condi- 
tions we choose a modified MS scheme, in which the heavy (anti-)quark is renormalized 
on-shell and the light quarks are renormalized using the standard MS prescription. The 
heavy quark masses are defined as pole masses. The subtraction for the renormalization of 
the strong coupling constant explicitly removes the heavy quark loop contribution to the 
gluon self-energy shown in Fig. 3 (i), see also ( B7|) in App. B. This leads to a fixed flavor 
scheme with the produced heavy flavor on the one hand and nif light flavors active in the 
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running of a s and in the parton evolution on the other hand [|T8j, [33| . The renormalization 
constants needed for the construction of the counterterms are then calculated to be 



Z m — 1 
Z 2 -l 

Z3-1 

Z g -1 

Zip 



9s» £ 

167T 2 



C F 3 



2 



(heavy quark mass), 



g 2 u~ £ 2 

" -Cp— (heavy quark field), 



16tt 

16tt 2 

9sV~ £ 



Sr.. 

2 2 2 
(2CU -&)- + -- 



327T 2 

Z„Z 2 Z 



.2 2 2 
Pot - 1: — 

1/2 



(gluon field), 

(coupling constant), 
(quark gluon vertex), 



with the QCD beta function (3q = (IICa — 2nj/)/3 for the riif active light flavors, and 
we have used the definitions 2/e = 2/e + 7^ — ln(47r) and 2/i m = 2m e /e with the Euler 
constant 7^. Z\p is determined to 0(g 2 s ) using the shown Slavnov- Taylor identity, and the 
quark photon vertex renormalization constant can be obtained from Z\f by either setting 



Ca = or by using the QED Ward-Takahashi identity Z?p D = Z 



The coupling 



constant g s in ([lSf) and in the matrix element calculations is the renormalized one in the 
counterterm method. We then of course use the renormalization group (RG) improved 
running coupling g s {y%) in the calculation, where fi r is the renormalization scale at which 
the subtractions are performed. We have checked that the procedure outlined above is 
completely equivalent to the one used in JTj| in the corresponding unpolarized calculation. 
In JTj| the renormalization is performed by replacing the bare mass and strong coupling 
constant in the Born cross section with the corresponding renormalized quantities (see 
IB] for more details). Their relation between the bare and the renormalized mass can be 



obtained by expanding mbare = Z m m in e using the Z m in flilf) . Analogously, the series in 
e of g!? arc = Z g g s using the Z g of (li"8|) in combination with the RG running leads to their 
relation for the strong coupling 



9, 



bare 



\ + 9 2 s (tf) 



m 

t4 



(19) 



327r 2 ^ \e 

where /3q is defined as po above, but with n/ = riif + 1 instead of riif flavors. The first 
term in the curly brackets corresponds to the usual MS prescription and the second one 
removes the heavy quark contribution to the gluon self-energy, as already mentioned. 
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(a) 



^RRnnr 1 — — mmmp 

(b) (c) 





(d) 



(e) 



(f) 



Figure 4: Feynman diagrams for the NLO gluon bremsstrahlung process *yg — >• QQg. Revers- 
ing the heavy quark lines yields the remaining graphs. In the unpolarized calculation the ghost 
contributions (e) and (f) have to be subtracted to cancel unphysical polarization contributions, 
see the discussion concerning Fig. |I|. 

The IR singularities of the virtual cross section are canceled by the soft part of the 
gluon bremsstrahlung contributions. The corresponding Feynman diagrams are shown in 
Fig. f| and the momenta assignment is 



7(fci) + g(k 2 ) -> Q{ Pl ) + Q(p 2 ) + g(k 3 



(20) 



In the calculation of this process we keep the kinematical invariants as defined in flT2| ) 
and introduce seven additional ones 



S3 = (fa +p 2 ) 2 - m 2 , s 4 = (fa + P1) 2 - m 2 , s 5 = (p 1 +p 2 )' 



-U5 



t> = (k 2 - fa) 2 
u& = (k 2 -pif 



m 



u' = (fa - faf 
= (fa - pi) 



(21) 



m 



Of course, only five of these invariants are independent for a 2 — > 3 process and hence 
there are many useful relations among the quantities in ( |i~2f) and (21), like 



X)! = — S — Ui — U-j , t' = — S — t\ — Uq , M 5 = t\ + U\ + S3 , 
S3 = S + M 6 + U 7 , S 4 = S + t\ + Mi . 



(22) 
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The real gluon bremsstrahlung cross section (R) can then be written as 

R 

F F G F [ dn £ \M R \ 2 , 




(23) 



r(l + e/2) 



,1+e 



27r(47r) 2+£ / 2 r(l + e) (s 4 + m 2 ) 1 + £ / 2 
d0 1 siii 1+e 1 / sur ft, 



'0 </o 

where we have used the standard 2-^3 phase space (see App. A). For simplicity we have 
already replaced the phase space integration over the (n — 4)-dimensional hat-space by "1" 
in (p3|) , since the hat parts do not give any contribution to the final result as discussed in 
Sec. 2 and App. A. Hence all hat terms in the squared matrix element |M^| 2 below have 
to be dropped before |Mr| 2 is inserted into (E3). G £ parametrizes the differences between 



the 2 — > 2 prefactor F £ defined in (15) and the 2 — > 3 prefactor, which is convenient for 
further considerations. The squared matrix element can again be split into an abelian 
and a non-abelian part 



\M R \ 2 = M^M R = E 2 £ g A s e 2 e 2 Q \2C F R QED + C A R OK 



(24) 



Concerning the angular integration J dQ e in (pj| ) we note that the Mandelstam vari- 
ables of (|l^) and ( plD are of three distinct types: Using the "set I" parametrization of 
the parton momenta in terms of the angles 9\ and 82 in App. A, the Mandelstam vari- 
ables t' and uq are of the [afe]-type a + fecos^i, whereas S3, S5, u' and u-j are of the 
[ABC]-type A + Bcos6\ + Csin^icos^- Here a, b, A, B, and C are functions of the 
angular-independent ([-]-type) Mandelstam variables s, ti, u\, S4 and of the heavy quark 
mass m. Extensive partial fractioning using the relations given in (|22|) reduces all phase 
space integrals to the standard form 

j(k,i) = j dQ e ( a + b cos 9 x y k {A + B cos fti + C sin Q x cos 6 2 )~ l , (25) 

which can be performed analytically. This reduction procedure can be completely autom- 
atized using a general set of transformation rules based on relations like in (|22"D , that is 
[ABC] = [■] + [ab] + [ABC]. As an example we demonstrate here two reduction steps for 
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1/ (s 3 u 6 u 7 ) using s 3 = s + u 6 + u 7 : 



1 f 1 1 1 \ , 1 1 1 

and = — , (26) 



s 3 u e u 7 s \u e u 7 u 6 s 3 s 3 u 7 j s 3 u 7 u 6 u 7 u 6 s 3 

with u'q = s + u e . In this way l/(s 3 u 6 u 7 ) can be completely reduced to the [a&][ABC] form 
required in (p5|), but the auxiliary [ab] variable u' 6 has to be introduced in our example. 
This always happens for two [ABC] variables in the denominator. Though the integrals 
with this kind of auxiliary variables are as straightforward as the others, they can also be 
easily avoided by introducing a different parametrization for the parton momenta ("set 
II" in App. A). In this particular set, u' and u 7 are of the [afe]-type instead of t' and u§ for 
"set I". So in our example above, l/(s 3 u 7 ) would be already of the [aft] [ABC] type using 
"set II" and no further reduction would be necessary. As a check for the correctness of our 
phase space calculations, we have proceeded in both ways. It should be remarked that a 
third conceivable parametrization of the momenta is not required in our calculations. 

The needed integrals j(*>0 are conveniently collected in JJT]. We have recalculated 
them and also the few additional ones that occur in the reduction method using solely 
"set I". The integrals are straightforward to calculate by using two "tricks": Firstly, 
most integrals can be derived from a basic one by partial differentiation with respect to 
the parameters a and A in (ESj). Note that one has to be careful in case of "collinear" 
parameters a 2 = b 2 or A 2 = B 2 + C 2 , though. Secondly, the transformation 

j(*,0 = j(W (a «_> A b _> + ~ hB ,C-> ~ bC ) (27) 

£ V y/WTC 2 VB 2 + C 2 ' K ' 

is often helpful. In particular one can transform 'A 2 = B 2 + C 2 " collinearities into 
"a 2 = 6 2 " ones in this way |35| . One can also prove Ie k,l \u§,u 7 ) = Ie' k \us,u 7 ,ti <-> Mi) 



using this transformation. 

In order to isolate the divergencies appearing in the soft S4 — > limit, which cancel 
the IR singularities of the virtual cross section, we examine the bremsstrahlung result in 



two regions: that for hard (S4 > A) and that for soft (S4 < A) gluons |18 |. Here the 
auxiliary quantity A is chosen small enough to be negligible in comparison to s, t\, u\ 
and m 2 . In the hard (H) region A effectively cuts off the IR singularities, so that only 
the M singularities remain. Thus no double poles will be encountered and one needs G e 
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in @3|) only to 0{e): 

Or = Or £ 

s 4 



167T 



1 + - [lE ~ ln(4vr) + In 



H {S4 + m 2 



:28) 



where we have also absorbed some additional factors into Gf for convenience. The hard 
gluon cross section then becomes 



a <J g-i,OK 

dt\dui 



H 



2 ~ 

CaF £ G £ E £ a s ae Q -H ok, pole 



H, 



OK,pole 



(p-l + l/p)- 

S 4 



hi 

Mi p 



Ui 



+ 0(1) 
Am 2 s 



P 



AH, 



OK ,polo 



2p- 3 + 2/p (t x 1 u x 



2m 2 s 



m 2 s 



(29) 



t X Ui 



where only the collinear pole contribution of the non-abelian OK part is shown. The hard 
abelian QED part is completely finite. The parameter p = 1 — s^/ui becomes one in the 
soft limit S4 — > 0, and one can clearly see the approach to an IR singularity proportional 
to the Born Bqed of (0). The finite contributions are too long to be presented here in 
an analytical form, but they can be found in our computer program, which is available 
upon request. Our unpolarized results agree with those of JT8|. 



Turning now to the soft gluon region, we find that one can write the Mandelstam 
variables in (^) in the soft limit s 4 — > as 

S3 = , t' = s 4 f , u' = s 4 u!_ , 

[60) 

U5 — —S + S4M5 , U 6 = Ui + > U 7 = ti + S4M7 , 

where the underlined quantities are finite dimensionless functions of the 2-^2 Mandel- 
stam variables (|T2|) and m 2 . In terms of these variables one can easily single out the IR 
singularities by collecting different powers in s 4 . Applying the s 4 — ► limit on the factor 
G e in (^), one obtains a s\ +£ dependence. Thus one only has to keep track of those parts 
of the squared matrix element exhibiting a 1/ s\ pole in the soft limit, since all other terms 
vanish for S4 —>■ 0. In this way one can easily derive the soft limit of Rqed and Rqk in 



S, 



QED 
SoK 



m 




QED, 









[I 


{-- 









QED, 



(31) 
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using the polarized and unpolarized Bqed and ABqed in QI^D , respectively, in agreement 
with PI . 



We rewrite the IR-divergent s 4 1+e -dependence in terms of the A-distribution[]: 



ds 4 f(s 4 )[g(s 4 )} A = ds 4 [f(s 4 )-f(0)}g(s 4 ), (32) 
'0 Jo 

where [(?(s4)]a is singular and f(s 4 ) is finite for s 4 — > 0. In particular we need the following 
identity 



-l+e 



A £ 

— 5(s 4 ) + [s 4 JA 



-1+el 



(33) 



This yields for the S4-integration of a function TL(s 4 ) with a soft pole s 4 1+6 S(s 4 ) and a 
finite T{s 4 ) part: 



ds 4 H{s 4 ) = / ds4 [s 4 1+£ S(s 4 ) + J r (s 4 
'0 



dsd 



A £ 



6(s 4 )S(s 4 ) + {S(s 4 ) - S(0)}s] 



-l+e 



^(34) 



<is 4 7i(,s 4 ) 



^-5(s 4 )S(s 4 ) + 6(s 4 - A)H(s 4 ) 



(34) 



where is the Heavyside step function. In the last step we have explicitly used that A is 
small enough to be negligible with respect to the 2^2 Mandelstam variables (|l2|) and 
m 2 , taking 5(A) ~ 5(0) ^ and J 7 (A) ~ ^(0) = as the finite limits. 



We can thus write the total gluon bremsstrahlung result by multiplying the hard cross 
section (|29|) with ©(54 — A) and adding the soft cross section obtained from (|31| ) with the 



-l+e 



replaced by 5(s 4 )A £ /e according to (|34). Using this replacement and performing 



the angular integrations dfl e the soft (S) gluon cross section is then given by 

s 



d 2 a (1) 

u °gj,OK 
dt\du\ 



d 2 d {1) 
a °g-y,QED 



n t? nS 772 2^ 2 Bqed 

Ca^s^e ^e a s ae Q n 



4 2 ti Mi 1 2 tii 
— + - In h In xln h - In 2 h 



Ml 



^ln 2 x + Li 2 (l- ' 



Ml 



-Li 2 1-— ] - 

XMi / V ^1 



2m 2 



8/3 



dt\du 



2C F F e G s £ E 2 £ a 2 s ae 2 ^^- 



2 1 2m 2 - S/i a . ..' 
-~ + l + ^-(l-S(x)) 



(35) 



3 If v = 1 + <i/s and u? = —Ui/s + t\ are introduced, then s 4 — > poles show up for u; — ► 1, i.e., 
l/s4 — > 1/(1 — u>). The singular u> — > 1 behaviour can then be treated with the usual +-distribution 
1/(1-«J)+ PI- 
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with 



S{k) 



— In x + Li 2 (x 2 ) - In 2 x + 2 In xln(l - x 2 ) - C(2) 



2(47r) 4 m 2 e 



Ge ^ l 67re£ (7 B -ln(4 7 r))/2 A _ ^( 2 ) £ 2 ^ ( 



A 2 



\i 2 m 2 



e/2 



S{s 4 ), (36) 



and where we have used f3 = a/1 — 4m 2 / s, x = (1 — + the dilogarithm function 
Li 2 as defined in |[36| , and the Riemann zeta function £(2) = 7r 2 /6. Adding the soft cross 
section fl35|) to the renormalized virtual cross section obtained from ([IB]) and flilf) removes 
all IR singularities in the latter, including the 1/e 2 poles. An additional 1/e pole in the 
OK part of (ftSQ will be eventually canceled upon adding the soft 5(1 — x) contribution 
of the mass factorization cross section (|37| ) discussed below. Our unpolarized results are 
again identical to those of |18[ . In addition we have checked that the abelian QED part of 
the polarized (and unpolarized) result is in complete analytical agreement with the NLO 
expressions for 77 — > QQ presented in [BO]. 

To obtain the final result for the gluon cross section, the remaining M singularities in 
the hard gluon cross section have to be removed as well. This can be achieved by the 
standard mass factorization procedure. To 0(ot 2 s ot) the reduced finite gluon cross section 
is given by |i~8| 



d 2 d {l) 
dtidui f 



d 2 d (l) 

(u 2 ) — 

dt\du\ 



(37) 



2tt 



dx 



P gg (x)-+Fgg(X } H 2 f ,fJ 2 



x 



2^(0) 



d 2 a 



91 



dtidui 



s 
h 



xs 



where fif denotes the factorization scale at which the subtraction is performed, and F 
represents the usual freedom in choosing a factorization prescription. In the MS scheme, 
which we use, F is given by 



Fij(x, (ij, /i 2 ) = Pij(x) ( lE - 1h(4tt) + In ^ 



(3? 



The Pij in ( |3"T| ) and ( 55 ) are the usual unpolarized (Pij) and polarized (APy) LO Altarelli- 
Parisi splitting functions |J7| 

1 1 

1 — x x 



9(1 - x - 8)2C A 



+ --2 + x(l-x))+P g 5 g (x) 
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AP gg (x) 



p5 M 



0(1 - x - 5)2C A 
/3q 



1 



5(1 



x 



1 — X 

+ 2C A hi6 



2x 



(39) 



Since we have regularized all soft singularities in our calculation by a small parameter 
A as outlined above, we have to stick to the same framework here to deal with the soft 
x — > 1 divergency of P gg consistently and cannot simply use the usual H — prescription 
1/(1 — x)+. In ( |59"D we have thus introduced another small auxiliary quantity 5 |18|, |38| . 
Of course, A introduced above and 5 are not independent, since they are related via the 
mass factorization by 5 = A/(s + £i). fio in ([39|) includes only the n\f light flavors as in 



Inserting (|5D|) in (R^) one gets schematically 

/ dx[5(l — x)A + 9(1 — x — 5)B(x)]x5(x(s + ti) + u{)C(xs, xt±,ui) 
Jo 



(40) 



5(s 4 )AC(s, t h mi) + 6(s 4 - A) 



-U\B> 



-Ml 



C 



-U\S —U\tx 



Mi 



S + tiJ \S + ti S + t\ 

where we have explicitly used the relation between 5 and A. Thus the contribution from 
mass factorization naturally splits into a soft and a hard part, which can be added to 
the corresponding cross sections. As already mentioned, the 2/e pole in the soft 5(1 — x) 
part of (|37]) cancels the remaining pole in the OK part of the soft cross section fl35|). In 
the hard cross section (|29|) the 2/e collinear pole is removed and one is left over with the 
finite contributions from the pole part obtained from the e-expansion of the prefactor. 
Examining this factor (|28|) and our choice for the MS factorization scheme in (0), it is 
easy to see that the reduced hard gluon cross section can be simply obtained from (|29|) 
by setting 



F £ G^ — > — 
s 2 



1 



and 



2 

£ 



In 



In 4 , (41) 



m 2 (s4 + m 2 ) 

where one can now perform the e — > limit. Again we agree with the unpolarized reduced 



hard gluon cross section of [y| . 



To complete the calculation of the gluon cross section, we now add the S(s^) mass 
factorization contribution in fl3"7D , see Eq. (flCf), to the renormalized virtual plus soft part 
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(V + S). We write the result in three parts using the usual abelian and non-abelian split 
and, in addition, separating off the part proportional to /3q/2. The latter piece vanishes 
if one identifies the renormalization scale with the factorization scale, i.e., [L r = ///, as is 
usually done: 

W\ _ a s^Q L r (J fA \ , r (f ,fA\,A)? 

dMsr J ~ r ' Q£D q ed ) + ^ (>* + l h + y L « 

The L A explicitly depend on the auxiliary quantity A. The polarized AL are presented 
in App. C and the unpolarized L are in complete agreement with those obtainable from 
App. A of |18|] and App. D of fl31| . The numerical treatment of the L A terms is discussed 
in App. C. 



• (42) 



To conclude this section we note that the presented results have been calculated for a 
detected heavy antiquark in the final state, because the heavy quark was integrated out 
in the calculations. But since all gluon matrix elements are symmetric with respect to 
Pi P2, the same gluon cross section holds for a detected heavy quark as well. On the 
other hand there is an asymmetry in the non-abelian part of the gluon cross section with 
respect to k\ <-> fe, since the outgoing gluon with momentum ks can only "couple" to the 
incoming gluon with momentum fc 2 , but not to the photon with momentum k\. 



5 NLO Light Quark Contribution 

In NLO one encounters a new type of subprocess with a light (anti)quark in the initial 

state 

f(*i) + q{k 2 ) - Q( Pl ) + Q(p 2 ) + q(k 3 ) , (43) 

which can be calculated along the same lines as the gluon bremsstrahlung contribution 
in the previous section using Eqs. (|2l| ) - (|23|) . The squared matrix element calculated from 
the graphs shown in Fig. |5| can be decomposed according to whether the photon couples 
to the heavy quark with charge cq (in units of e) in the "Bethe-Heitler-graphs" (a) and 
(b), or to the light quark with charge e q , as for the "Compton-graphs" (c) and (d): 

\M q \ 2 = A^ = E e9 y^[(? Q A 1 + e 2 q A 2 + e q e Q A 3 \ , (44) 
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(a) 



(b) 





(c) 



(d) 



Figure 5: Feynman diagrams for the NLO light quark initiated process — > QQq. 

where A% denotes the interference contribution of both types of processes. Notice that 
since we have now only one boson in the initial state, the photon, only one factor E e 



appears in (0). Since this production mechanism appears for the first time in NLO 
and there are no gluons in the final state, we do not encounter IR singularities in the 
calculation. All single poles can be solely attributed to collinear configurations, and thus 
an expansion to 0(e) is sufficient here as for the hard gluon cross section in Eqs. fl2~8|) and 



The phase space integrations and the preceding partial fractioning proceed just as was 
explained in Sec. 4 and so we can immediately quote the results here 



dt\dui 

^pole 
AA polc 
^pole 

AA p ° lc 



^F £ G»E £ a 2 s ae 2 -Af c + 0(1) 



e 



2(p-l) + l/p 



VP 



hi 

t\ 1 U\ 

— + —P 

ui p h 



r-2 + 2/r 



U\ Am 2 s 

ti tiU\ 

2m 2 s 



2m 2 



T 



tH{ + u{ 2m 2 



m 2 s 



(45) 



where again only the collinear pole contributions are given and e 2 denotes Cq, e 2 , and tQe q 
for i = 1, 2, and 3, respectively. The interference contribution A 3 is completely finite. The 
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new parameter r = 1 — 8±/t\ becomes one in the limit S4 — > 0, exactly like p = 1 — S4/U1 
already introduced in (fZ5|). A^ ' 6 stems from 1/t' terms in the matrix element (0). £' 
becomes zero for the collinear configuration k 3 = /c 2 (l — x). In "set I" this zero shows up 
as a = — b with a = puis^/[2(s^ + m 2 )]. Note that A± olc becomes proportional to Bqed 
in the limit — > 0, and thus S4 — ► 0, for x — > 1. Similarly, 74 2 ole originates from 1/V 
terms, w' becomes zero for £3 = &i(l — x) and in "set I" one finds A = — \/B 2 + C 2 with 
A = Tt 1 s 4c /[2(s 4 + m 2 )}. In the limit x -»• 1, i^ olc becomes proportional to Aqed, which 
is defined in (|49|) below. The complete expressions for this subprocess are too long to 
be given here but can be found in our computer program. Our unpolarized results fully 



agree with those of [18 



Again the mass factorization procedure removes the collinear singularities. To 0(a 2 a) 
the reduced finite quark cross section is given by fl8| 



2aW 



crer, 



dt\du\ 



d 2 a {1) 

7 (/j 2 ) 
dtidui 



(46) 



a* 



/ dxi 

2vr ' 



/ c/x 2 



+ F gq (xi, fl 2 f , fl 2 

£ 



P q j{X 2 ) - + Fg 7 (X 2 ,/ij,/i 2 ) 



Xl 



x 2 



'd 2 ~afj' 






dtidui 


G:: 








-> x 2 s 


dt\du\ 







and for light antiquarks the analogous relation with q — > g holds. The first subtraction 
in (|46| ) corresponds to the collinear configuration in the Bethe-Heitler part, whereas the 
second piece refers to the collinear Compton contribution. The quark-gluon 

l + (l-x) 2 " 



Pgq - Cp 



X 



AP gq = C F (2-x) 



(47) 



and the photon-quark splitting functions 

P qi = e 2 q N c [x 2 + (1 - a;) 2 ] 



AP 97 = e 2 N c {2x - 1) 



(48) 



in (f46|) can be obtained from ||37|| . The corresponding antiquark splitting functions are 
identical and Fij was already specified in (|38D in the MS scheme. 



It should be noted that the subtraction term proportional to P qi in ( f46|) implicitly 
introduces the quark content of the real (on-shell) photon. Clearly, one cannot obtain 
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a factorization scheme independent result taking into account only the "direct" point- 
like photon contribution without adding the corresponding "resolved" cross section which 
probes the parton content of the photon. This is made evident by the appearance of 
F qi in fl46|) , allowing for arbitrary redefinitions of the factorization scheme (i.e., of the 
photonic quark densities in NLO), which can only be compensated by the NLO resolved 
contributions. Since the spin-dependent resolved cross section has not been calculated 
in NLO yet, it has to be estimated in LO. A further complication arises here, because 
the parton content of longitudinally, i.e., circularly, polarized photons is experimentally 
completely unknown for the time being, and one has to rely on realistic models [ |3TS[ 
when estimating the size of the resolved contribution. However, it has been demonstrated 
in [[RJ that even for large spin-dependent photonic densities, the "background" from 
resolved photon reactions should be very small for all experimentally relevant total or 
differential cross sections. In particular this is the case at fixed target energies, as for 
COMPASS. Only for the total charm production spin asymmetry at collider energies the 
resolved contribution can be as large as 30% fl6fl. But in this kinematical region, charm 



production anyway suffers from large statistical errors and appears to be unmeasurable 



at the polarized HERA option [|12 



It should also be remarked that for NLO photonic parton densities, unpolarized [[|(J 



as well as polarized [[|l|] ones, often the so-called DIS 7 factorization scheme rather 



than the MS prescription is used, since it provides a better perturbative stability between 
LO and NLO quark densities. In this case one either has to transform the the densities 
back to the MS scheme p0|, [4T| before using them in the calculation of the NLO resolved 



contribution or one has to use the appropriate DIS 7 expression for F„ 7 in (Hq), see the 



Appendix of [42 



To calculate the factorization contribution for the Compton part proportional to e 2 q in 
(fPJ), the Born cross section for the qq — > QQ process in n dimensions is required. Only 
one diagram contributes here, see Fig. |6|, which can be easily calculated using the general 
2^2 phase space expression of (|T5|) and the result reads 

d 2 a {0) r 

- f ^W c 9sAqed ' 
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Figure 6: The LO quark-antiquark annihilation process qq — > QQ. 



t? + u? 2m 2 e 

H h - 

s 2 s 2 

t 2 + m 2 2m 2 £ 
i 2 ~ + 2 



A QED = + ^ + i , (49) 

3 



The mass factorization in ( f46|) is again performed with the MS prescription (p8[), and 
the finite reduced quark cross section can be obtained by applying Eq. ([0]) to (fS|). Our 
results fully agree with |18[ in the unpolarized case again. 



Finally, it is important to point out that A\ and A 2 in stay unchanged for pi <-> p 2) 
whereas A3 changes sign. Thus if one wants to observe a heavy quark instead of a heavy 
antiquark, one can use [cqAi + e 2 A 2 — e q eQA 3 ] in fl44]) with the same expressions for the 
Aj. In addition charge conjugation implies that 

da (<77 — ► Qj = do 1 (57 — ► Q) and da (777 — > Q) = da (57 — ► Q) . (50) 

Thus one can use the same Aj for the contribution due to a incoming antiquark in (|43|) 
as well, taking into account a negative sign for A 3 . Note that the sign change of A 3 also 
implies that A 3 does not contribute to the total cross section ([51]) below, since the result 
cannot depend on whether the heavy quark or heavy antiquark is integrated first. 

6 Hadronic Cross Sections and Numerical Results 

Let us now turn to some phenomenological aspects. The total photon-parton cross section 
can be expressed in terms of scaling functions in both the unpolarized and polarized case 
(i = g, q, q) MM- 

^(s,™ 2 ,/^/,/^) = / d{-t x ) [ _ d(-u l ) d ai ^^ U ^ ( 51 ) 



aa 



m 2 



ffl ( V ) + 4tt« J /« (rj) + ( V ) In & 
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where rj = s/ (Am 2 ) — 1 and (3 is defined below Eq. (|36|) . fff and fj^\ f { stand for the LO 
and NLO corrections, respectively^ This coefficient functions can be further decomposed 
depending on the electric charge of the heavy and light quarks, eg and e q , respectively: 

ftniv) = e QC g7 {v) , (52) 
frriv) = CQCqiiv) + e%i(v) , (53) 

with corresponding expressions for the f iry . Note that the interference contribution pro- 
portional to tQt q drops out in f qi as discussed in Sec. 5. 



The behaviour of the spin-dependent coefficient functions ( p2j ) and fl53|) has been 
already shown and discussed in detail in for the conventional choice \x$ = fi r . Here 



we just want to point out that for fif ^ fi r one receives an extra contribution from L^p 
in (^2[), see also Eq. ( |C5| ) in App. C. This can be easily accounted for by adding an 



appropriate term to the NLO gluonic coefficient function 

/«W=45«W+^'Wln| , (54) 

and, of course, by using \i r as the scale for a s . Notice that this is equivalent to the 
replacement 

ot s (i4) -> a s (i4) ( 1 + a s (/4)^- In ) , (55) 



in Eq. ([5TD, keeping only terms up to O(a 2 s o). We will study the effect of varying /i/ and 
/i r independently on the total hadronic heavy flavor photoproduction cross section given 
by ^ 

a® p (S,m 2 ,fi f ,fi r ) = ^2 / dxcr(xS,m 2 ,fi f ,fi r )f(x,fi 2 ) (56) 

in detail below. S in ([56]) denotes the available photon-hadron c.m.s. energy and the / 
are the (un)polarized parton distributions. 

Before that let us first of all recall the relevant formulae for calculating differential 
single- inclusive heavy (anti) quark distributions. We denote the momenta in the photon- 

4 Note that we still use the "tilde notation" as a shorthand to denote both the longitudinally polarized 
and unpolarized cross sections simultaneously. 
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hadron cross section by 

7 (fci) + H(K 2 ) - Q(p 2 ) [Q( Pl )\ + X, (57) 

and use the following hadronic invariants for the observed heavy antiquark 

S = (k 1 + K 2 ) 2 = ~ , T x = (K 2 -p 2 ) 2 -m 2 = t ±,U x = (k.-p^-m 2 = u x , (58) 

where we have introduced the momentum fraction x in k 2 = xK 2 to relate the hadronic 
to the partonic variables in Eq. flTjp. For an observed heavy quark one would exchange 
Pi *-* P2 in Q58|) and there would be a t x ^ u x crossing. 

The hadronic single-inclusive heavy (anti) quark cross section reads 



dTidUi ^—^_Jt . J y ,r ~ f/ dt x du x 

g,g,q " x mm 

and the lower limit x m i n of the integration is determined from 



dxxf(x^j)—^- , (59) 



s 4 = xS + xTx + U x > A => x min = — ^ . (60) 



For the actual integrations it is convenient to change the variable from x to in ( |59| ) 
with the limits A < s 4 < s^ ax = S + T x + U x : 

/ dx= s + T / dSi ' ^ 
where A is the cutoff introduced above to separate the hard and the soft gluon cross 
sections. The soft plus virtual gluon cross section proportional to 5(s4) has to be evaluated 
with elastic kinematics (s + t x + u x = 0). However, for numerical purposes we rewrite 
the In* A/m 2 (i = 0, 1, 2) terms in da s+v into integrations over s 4 as outlined in App. C. 
In this way the soft plus virtual and the hard parts of the gluonic cross section can be 
directly added. One can always set A = for the light-quark induced subprocess. 



The differential heavy (anti) quark cross section ( |5T?[ ) should be expressed in variables 
more suited for experimental measurements: 

transverse momentum/mass Xt = ^ T , m\ = m 2 + p\ — ' " 



Pt ° 

xtanh— = - In — 
E 2 T 

PL l_ T\-Ui 

PT X Ps s 



rapidity : y = artanh^f = - In — ^ , (62) 

E 2 T 1 

Feynman 



26 



where @s = yl — Am 2 /S. The energy and the longitudinal momentum of the heavy 
antiquark are given by = m^coshy and pi = m^sinhy, respectively. j>r = \pt\ is the 
absolute size of the transverse momentum and p™ ax = p™ ax = \/~S(3s/2. y and of the 
observed Q in (|62|) are defined in the hadronic c.m.s. with the forward direction (y > 0) 
along the incoming photon, i.e., 

T x = -^fSm T e- y = -VSp^ix ~ x F ) , XJ X = -^Sm T e y = -VSp^ix + x F ) , (63) 



where % = \[xp + (mr / p™ ax ) 2 . 

The Jacobians needed to express ( |59|) in the variables ( |62|) are 

dTidUi = Sdm T dy = —dm^dxp , (64) 

and cfmy = 2x r (p™ ax ) 2 <i£'r, etc. By integrating the variables in ( |62]) over the appropriate 
limits 

/ ' 2 1 T / 1- AS / 4 cosh 2 » „ 

5 / dm T I dy = S / cfa/ / <im T , 

Jm 2 i-arcosh^S- ./ 1 l n *±|fi. J m 2 

S(a_o2 „2 • 



4 m 



s ^ = s f dXF ** (65) 



/s V " ; 



1-^ X 7-1 7m 2 X 



the total cross section (1561) is of course recovered. 



Finally it should be noted that experiments do not determine the (differential) longitu- 
dinally polarized cross section (d)Acx itself, but rather the corresponding spin asymmetry 

Q _ (d)Aa 

In (|66|)) which is nothing but the counting rate asymmetry for the two possible helicity 
alignments of the incoming photon and hadron in analogy to Eqs. @ and (|7|), the exper- 
imental normalization uncertainty and other systematical errors conveniently drop out. 
However, in the following we will concentrate on the polarized cross section itself, since 
we are mainly interested in the influence of the spin-dependent NLO corrections. The cal- 
culation of the spin asymmetry (|66| ) would introduce additional theoretical uncertainties 
associated with the unpolarized (differential) cross section. 
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Equipped with the necessary technical framework, we now turn to some numerical 
applications. Unless otherwise stated we use here the GRV |43| and GRSV || "standard" 



set of unpolarized and longitudinally polarized parton distributions, respectively, in our 
calculations. The strong dependence of the results on the chosen, experimentally poorly 
constrained polarized gluon distribution Ag (and to a lesser extent also on the unpolarized 
gluon g) has been already demonstrated in |20] for the case of the total charm production 



cross section and the corresponding spin asymmetry. Of course this sensitivity in turn 
implies that such a measurement is particularly suited to pin down Ag more precisely. 



Unfortunately we have no data so far, but in the near future COMPASS is going to 
measure the total charm spin asymmetry A c lp with sufficient accuracy [13|, |2(J . Therefore 
we mainly focus on the kinematical range accessible by COMPASS in our analyses below, 
i.e., y/S = \J S 1P = 10 GeV. It is currently under scrutiny whether it is physically 
feasible and sensible to run HERA in a polarized collider mode in the future |P2]| , and 



therefore we either show or comment on the corresponding results at HERA collider 
energies y/S ~ 100 — 300 GeV as well. 



In |20| we kept the renormalization and factorization scales equal at the common 
choice fif = fi r = 2m for the total charm spin asymmetry. Here we investigate the 
theoretical uncertainty induced by varying fif and fi r independently in the range /i 2 , fif. = 
m 2 , . . . ,4.5m 2 with m = 1.5 GeV for the charm quark mass. Fig. 7 shows the deviation 
R of the LO (a) and NLO (b) polarized total cross sections Aa^ p (]56|) from the results 
obtained for the choice /x 2 = /x 2 = 2.5m 2 , i.e., 

„ = A ^p(/4/4) - Aor 7p(^ = /4 = 2 - 5m2 ) (67] 

A^ p (/z 2 = /i 2 = 2W) 1 j 

for v^S* = 10 GeV in percent. The contour lines are in steps of 5%, and at the base of the 
plot the line corresponding to the usual choice fif = fi r is shown for convenience. As can 
be inferred from comparing Figs. 7(a) and (b), the scale dependence has been drastically 
reduced in NLO in the entire range for fif and fi r , which underlines the usefulness of the 
NLO results. Moreover, in NLO the choice fif — fi r is approximately on the contour for 
R = 0, and R is flattest for large fif and fi r . This a posteriori motivates and justifies our 



choice of scales, fif = fi r = 2m, in PD[ . For reasonable changes of fif and fi r in Fig. 7, 
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R (LO) R (NLO) 




Figure 7: Renormalization /i r and factorization /ij scale dependence of the LO (a) and NLO 
(b) R as defined in Eq. flo7|) in percent for \fS = 10 GeV. fi r and fif are in units of the 
charm quark mass m c = 1.5 GeV. The contour lines are in steps of 5% and for convenience 
a line corresponding to the usual choice fx r = [if is shown at the base of the plots. 

the polarized total charm production cross section (|56|) varies by about ±10% in NLO as 
compared to about —15% to 35% in LO. It should be noted that one finds very similar 
results also for a higher c.m.s. energy of, e.g., y/S= 100 GeV. 

Fig. 8 illustrates in a similar fashion as in Fig. 7 the dependence of the polarized 
and unpolarized total charm photoproduction cross sections fl56|) and the corresponding 
longitudinal spin asymmetry (|66|) on the mass of the charm quark for two values of \^S. 
We vary m around our standard choice m = 1.5 GeV by ±0.2 GeV, and R m in Fig. 8 
is defined in analogy to Eq. flBTp . As is expected, far above the production threshold for 
y/S = 300 GeV ^> 4m 2 , the dependence of the quantities in Fig. 8 on the precise value 
of m is strongly reduced as compared to the results obtained for \/S = 10 GeV. For 
a reliable extraction of Ag by COMPASS, the mass uncertainty is far more important 
than the scale dependence in Fig. 7, since the experimentally relevant spin asymmetry 
A c varies by as much as 30% in the shown mass range. A good determination of m 
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Figure 8: Mass dependence of the (un)polarized total charm photoproduction cross section 
(j56|) and the spin asymmetry (^) in LO and NLO in terms of the ratio R m as defined in the 
text for t/S = 10 GeV (denoted by a cross) and 300 GeV. 

is thus mandatory for a meaningful determination of Ag at low energies, not too far 
above threshold. In addition, we have already stressed in |[20|| , that for a determination 
of Ag at fixed target energies further complications arise also from our poor knowlegde of 
the unpolarized gluon distribution for x > 0.1. Unfortunately, a measurement of A c lp at 
collider energies, where the theoretical uncertainties are much better under control than 
at low energies, appears to be not feasible, since A^ p is at best of the same size as the 



expected statistical errors for such a measurement [16|. This also does not improve for 
Pt or y differential charm distributions. 

In Fig. 9 we turn to the longitudinal spin asymmetry A b lp for total bottom quark 
photoproduction (m = 4.5 GeV) in NLO for HERA energies for four different sets of 
polarized parton distributions || B], |TT| which mainly differ in Ag. The results obtained 
for the different sets of parton densities are well separated and sensitive to the different 
Ag, but A b is extremely small. Since A^ p already appears to be unmeasurable at HERA, 
the prospects for a meaningful measurement of A h seem to be not very promising at the 
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Figure 9: The longitudinal spin asymmetry A b for total bottom photoproduction in NLO 



in the HERA energy range for m = 4.5 GeV, /!/ = fi r = 2m and four different sets of 
NLO polarized parton densities || |], [I1J. The "error band" is an estimate for the expected 



statistical accuracy 5A for such a measurement at a polarized HERA (see text). 



first sight, since bottom cross sections are smaller due to the larger b quark mass and 
the smaller heavy quark charge (e{,/e c ) 2 = 1/4. However, b quarks are experimentally 
much easier to detect, e.g., through their longer lifetime (secondary vertex tag), which 
might compensate these shortcomings. The shaded band in Fig. 9 illustrates the expected 
statistical accuracy for such a measurement at HERA estimated via 

1 1 



6 A" ~ 

7P — 



(6* 



assuming a polarization P e p of the electron and proton beams of about 70%, an integrated 
luminosity of C = 500 pb -1 [O, and an optimal detection efficiency of e& = 0.05 E4 . 



Finally, let us turn to some results for differential distributions. Although their exper- 
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Figure 10: i T differential polarized anticharm photoproduction cross section dAa^ p /dxT in 
LO and NLO according to Eqs. (|9])-(||) for y/S = 10 and 200 GeV and two sets of polarized 
parton distributions [§, [11]]. The bands with forward (NLO) and backward (LO) slanted hatches 
correspond to the uncertainty due to independent variations of fij and ji r (see text). LO results 
are denoted by stars, and the "GRSV std." curves for y/S = 200 GeV are multiplied by 1/8. 

imental relevance seems to be remote, apart from p? and y acceptance cuts, a comparison 
of the LO and NLO distributions is of theoretical interest to understand in which kine- 
matical regions the corrections are most relevant. 

In Fig. 10 we show the rapidity integrated polarized cross section in LO and NLO as 
a function of X? according to Eqs. (|59|)-(|65D. Two values of \/S (10 and 200 GeV) and 



the GRSV "standard" || and DSS "set 3" |Tl[ polarized parton densities are used. For 
the GRSV results the theoretical uncertainty of varying /i^ and \j? r independently in the 
range a{p\ + m 2 ) with a = 1/4, . . . , 4 is illustrated by the bands with forward (NLO) and 
backward (LO) slanted hatches. All curves are calculated for the choice a = 1 and the 
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Figure 11: Rapidity y differential polarized anticharm distribution dAa^ p /dy at = 
lOGeV in NLO and LO according to Eqs. (||9])-(|65]) and using fi f = /i r = 2m. The light 
quark induced NLO contribution is shown separately. For comparison the dotted curve displays 
the NLO unpolarized distribution scaled down by a factor 8. 

LO results are marked by stars. The NLO corrections are sizable, but the NLO shape is 
very similar to the LO one. Note that the large corrections for the GRSV || densities are 
to a large extent due to the differences between the poorly constrained LO and NLO Ag. 
We have made a similar observation for the total cross section in [20]; using the GS |J 



densities (not shown) this effect would be even more pronounced. The GRSV and DSS 
curves only differ in size, due to the much smaller DSS "set 3" gluon, but not in shape. 
As for the total cross section in Fig. 7 it turns out that the scale uncertainties are reduced 
in NLO. For example, at xp = 0.1 the NLO result varies by —20% to 55% whereas the 
LO result varies by —35% to 75%. For y/S = 200 GeV the scale dependence is similar in 
NLO, but even worse in LO. 



In Figs. 11 and 12 we show the c.m.s. rapidity y and Xp differential polarized anticharm 
photoproduction cross sections according to Eqs. (|59|)-(|65|). Since pp is integrated over 
the entire kinematical range, we choose fif = fi r = 2m as the hard scales here. The 
distributions are asymmetric in y and xp and the heavy quark is dominantly produced 
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Figure 12: As in Fig. 11, but now as a function of xp- The NLO unpolarized result is divided 
by a factor 11 here. 

"backward" with respect to the incoming photon, i.e., in the direction of the proton. The 
NLO results are always larger than the LO ones and deviate in shape. In both figures 
the unpolarized distributions, scaled down to the size of the polarized one, are shown for 
comparison. The genuine NLO contribution with light quarks in the initial state is shown 
separately and appears to be negligible in the entire y and xf range. 



7 Summary 



To conclude, we have presented the details of the first complete NLO QCD calculation 
of heavy flavor photoproduction with longitudinally polarized beam and target. We have 
provided all relevant intermediate steps of our calculation, in particular we have given 
complete analytical results for the soft plus virtual gluon cross section. A compact nota- 
tion was introduced to present both the unpolarized and polarized results simultaneously, 
and whenever possible we have compared our results to the existing unpolarized calcula- 
tions and found complete agreement. Similarly, for the abelian part of our unpolarized 
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and polarized results, which can be compared to Refs. |3(J analytically. 

As phenomenological applications of our results we have first explored the theoretical 
uncertainties due to the independent variation of the factorization and renormalization 
scales and due the unknown precise value of the charm quark mass. It was found that the 
scale dependence is much reduced in NLO, which clearly demonstrates the usefulness of 
our NLO results in future determinations of the polarized gluon density Ag, for instance, 
by the COMPASS experiment. It was critically discussed that the value of the charm 
quark mass is one of the major uncertainties in a measurement of Ag at fixed target en- 
ergies. NLO estimates for the total bottom quark spin asymmetry accessible in a possible 
future polarized collider mode of HERA were presented. Finally, we have presented for 
the first time xt, y, and xf spin-dependent differential single anticharm photoproduction 
cross sections in NLO. Although their experimental relevance seems to be remote, our 
differential expressions are useful for pt and y acceptance cuts for upcoming "total" cross 
section measurements, in particular at COMPASS. 
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Appendix A: Kinematics and Hat Phase Space 



The phase space calculation for the 2^3 processes is conveniently performed in the 



c.m.s. frame of the two outgoing unobserved partons [24, 31] (we follow here as far as 



possible the notations in App. B of pl| ): 

ks = (u 3 , , u>s sin 9\ cos 9 2 , u 3 cos 9i , k 3 
pi = (Ei, pf, -u 3 sin 6i cos 9 2 ,-u 3 cos 9i,pi) . (Al) 

Since k 3 +pi = 0, only one (n — 4)-dimensional vector k = k 3 = —p± remains and similarly 
fcf = —Pi- The ^-components need not to be specified, since the matrix elements do not 
depend on them, and hence the x-integrations can be trivially performed. The other three 
observed momenta k±, k 2 , and p 2 can be oriented in such a way that they lie in the y — z 
plane with vanishing hat components. There are three sets depending on which vector 
is chosen to point in the z direction. As outlined in Sec. 4, the whole calculation can be 
performed by using only one of these sets. "Set I" is given by 

k\ = (wi, 0, |p| sin 1 ]/, \p\ cos* — u> 2 , 0) , 

k 2 = (u 2 , 0,0,w 2 , 6) , (A2) 
p 2 = (E 2 ,0, |pl sin^, |p1cos^,6) , 

where the kinematical quantities in ( |Al| ) and ( |A2| ) are determined from on-mass-shell 
constraints and momentum conservation fl3T||: 



S + U\ S + t\ S4 

U l = n , . « » ^2 = — ~ 1 == , UJ 3 



2\/s4 + m 2 ' 2V^4 + rn 2 ' 2^4 + rn 2 

5 4 + 2m 2 h + ui + 2m 2 

1 = O / i 2 ' h2 = o / | 2 ' ( A6 ) 

\/(ti + ui) 2 — Am 2 s tis 4 — s(ui + 2m 2 ) 

\P\= " r^ =^ » cos*- 



In a second parametrization of the momenta ("set II") u' and u-j in (^T|) are of the [ab]- 
type instead of t' and u 6 in set I. The combined use of these two sets, depending on which 
type of Mandelstam variables in the denominator has to be integrated, is sufficient to 
avoid any appearance of auxiliary quantities like u' 6 introduced in (|26|) . In set II we have 

ki = (uj 1 , 0,0, ui,6) , 
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k 2 = (<j 2 ,0, |plsintf, |p|costf-£Ji,6) , (A4) 
p 2 = (E 2 ,0, |pl sin^, |p1cos^,6) , 



with the same relations as in (|A3|) except for 

, T , UiSi- s{t l + 2m 2 ) 

cosW = : . (A5) 

(s + Mi) a/ (ii + u\) 2 — Am?s 

The derivation of the two-to-three body phase space formula dPS 3 is standard, and we 
concentrate here only on the new aspects due to the the additional (n — 4)-dimensional 
hat-space integration in the polarized case. The calculation of dPS3 is facilitated by 
introducing a pseudoparticle with momentum p = pi + ks, i.e., the sum of the momenta of 
the two unobserved partons. dPS3 can then be separated into a 2 — > 2 and a 1 — > 2 phase 
space. Only the latter "decay" of the pseudoparticle into unresolved partons depends on 
the hat-space with the choice of coordinates explained above. This non-trivial integration 
is then given by 

d n k 3 d n p 1 5 + {kl)5 + {p\ — m 2 )5^ n \p — k 3 —pi) . (A6) 



Using ( [AT] ) and (|A"3|), i.e., p = (\/s4 + m2 , 0, 0, 0, 0), and the fact that the matrix elements 
depend only on k 2 , we can evaluate (|A6|) easily by integrating k 3 , p\ and the angles of k 



71 6 4 / / in ,n • 2 



r(n/2-2)8(s 4 + m 2 ) 3 / 2 

^2(p)n/2-3 



JO 



.-I 







4(s 4 + m? 



d9id92sm #isin#2 
sin 2 9\ sin 2 82 — k 



-.n/2-2 n-3 pir pn pi n/2-3 

' ' dM0 2 sin™" 3 ^ sin n " 4 6 2 / dx- 



r(n/2-2)[4( S4 + m 2 )]«/ 2 -iio Jo Jo 

_ W 2 - 2 r(n/2 - 1) ,r 3 /■ 

" r(n-3) ( S4 + m 2 )W2-i / dU "- 4 2 ' 

with the following abbreviations for the remaining integrations 

/PIT P7V 
dVL n ^ =11 dd x dQ2 sin™" 3 #1 sin™" 4 6 2 , (A8) 
Jo Jo 

1 /-l ™n/2-3 

J = — ; / dx , . (A9) 

S(l/2,n/2-2)i v/T^ 



Furthermore we have used the definition 

x = k lk 2 max = \ — . y 2 . (A10) 



4(s 4 + m 2 )fc 2 
s 4 sin 2 6*i sin 2 # 2 
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Comparing ( |A7| ) with (BIO) in Ref. pif , we see that the dependence on the hat 
momenta has been completely absorbed into the additional integral X. Thus we can 
schematically write for the full 2 — > 3 phase space formula 

dps 3 (0i, e 2 ) = dPS 3 , unp (#i, e 2 ) x i , (Aii) 

where (one should keep in mind that (B14) in [RTT| already includes the flux factor, etc.) 



n i „ 

mn , n n s , , 1 (4vr)- n ft l u 1 -m 2 s\ 2 s^T 3 , , 



Notice that ([A~^) - flA~T2"D agree in the limit m 2 — > with the result presented in [p6 |. 



The vast majority of terms in the matrix elements do not depend on k 2 and the rest 
is proportional to k 2 , so that one finds only two cases (n = 4 + e): 

i-i dPS 3 {6 1 ,e 2 ) = aps 3tWp {e 1 ,e 2 ) , 

X-P =► dPS 3 (fli,fl 2 ) = dPS 3 , unp (g 1; g 2 ) e gl^i^ , ( A13) 

The IR poles in our calculation stem from terms diverging as 1/s 2 for — > 0. Such poles 
are canceled by the factor s 2 in (|A13|) . The M poles in the matrix elements have the 



collinear structure sin 1+e 1 /(l — cos#i), and again we find that due to the factor sin 2 9% 
in ( | A 1 3| ) one gets finite results. Thus all hat integrations turn out to be infrared and 



collinear safe. Due to the extra e in ( |A13| ) all hat contributions are then of 0(e) and drop 



out when the e — > limit is taken in the end. Note that the heavy quark mass plays the 



role of an infrared regulator here. Collecting the S4-factors in ( |A13[ ) and ( |2"3"D one gets 
s^ +£ /(s4 + m 2 ) 2+£ l 2 . For m — > one then has s]^ £ ^ 2 } which is not sufficient to cancel the 
I/S4 contributions. Thus in the massless case one picks up extra finite contributions from 
the hat-space integration with infrared poles 26|, whereas collinear safety still holds. 



Appendix B: Tensor Integrals - Some Remarks 



When performing the Passarino-Veltman decomposition |32[ various scalar integrals ap- 
pear which at a first glance are not listed in [|3T] . However, they can always be cast in the 
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form of | 3"Ifl . For example, shifting the loop momentum q — > q + p 2 immediately yields 
C (pi,-fci - k 2 ,0,m,m) = C (p 2 ,pi,m, 0,m) . (Bl) 



In less simple cases one can easily find the necessary relations by inspecting the Feynman 
parametrization of the integral. For the three point functions Co(q±, q 2 , mi, 777.2, m^) the 
Feynman parameter integral has a denominator [q 2 — C] 3 with C = —abq\ — ac(q± + q 2 ) 2 — 
bcq\ + am\ + bm\ + cm 2 , where a, b, c are functions of the two Feynman parameters with 
a + b + c = 1. We can then prove, for instance, 

C (pi - fci, -fc 2 ,0, m,m) = C (pi, -fci,0,m, m) . (B2) 

by simply inserting the momenta and masses in C and interchanging b <-> c. It should be 
noted that exploiting the freedom to re-assign the parameters is also essential for explicitly 
calculating the set of basic scalar integrals. 

When decomposing the tensorial four point functions for the box graphs in Figs. 3 (a) 
and (b), one has to keep the rather lengthy intermediate expressions as short as possible. 
For the calculation of the QED-like box Fig. 3 (a), one can show that eight of the 
twenty-three scalar coefficients are not independent: 

D u = D l2 + D 13 , D 2A = (D 21 + D 22 - D 23 )/2 , 

D 25 = (£> 2 i - £22 + D23V2 , D 31 = -2D 32 - 2D 33 + 3^36 + 3D 37 , 

D u = ~D 32 - D 33 + 2D 36 + D 37 , D 35 = -D 32 - D 33 + D 36 + 2D 37 , (B3) 

Ano = (~D 32 - D 33 + ^ + D 37 + £>3 8 + D 39 )/2 , Dan = D 312 + D 313 . 

To prove this, one starts by replacing q — > — q in (^) and with — g 4 = q± + q 2 + q 3 finds 

D (q 1 ,q 2 ,q 3 ,m 1 ,m 2 ,m 3 ,m A ) = D (q 4 , q 3 , q 2 ,m 1 ,m i ,m 3 ,m 2 ) . (B4) 

We will abbreviate this result as Dq = Dq. Each additional power of the loop momentum 
in the numerator introduces an extra factor (—1), so that 
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By writing down for example the vector decomposition and comparing the coefficients, 
one obtains 

D u = D u , D 12 = D n - D 13 , D 13 = D n - D 12 . (B6) 

In the case of the QED-like box one encounters D^ = £)y(pi, — ki, — k 2 , 0, m, m, m) 
and = Dij(p 2 , —k 2 , —k%, 0, m, m, m). But the QED-like box is symmetric with respect 
to k\ <-> k 2 and p\ <-> p 2 . Thus D^ = and one obtains the first relation of ( |B"5| ) 
from ( |B6| ) and, analogously, the other relations are derived. For the non-abelian box 
in Fig. 3 (b) on the other hand, one has Dij = Dij(— ki,p±, — k 2 , 0, 0, m, m) and D^ = 
Dij{p2i — k 2 ,Pi, 0, m, m, 0) and furthermore the ki <-> k 2 symmetry is lost as discussed 
at the end of Sec. 4. Thus here one finds no simple relations between the D^ and the 
Dij that could be exploited to reduce the number of coefficients. One can, of course, 
check these general arguments explicitly, e.g., for arbitrary momenta and masses one 
finds D12 —Dyi — D%3 — D13, see ( B6|) . But this difference only vanishes in the QED case, 
whereas in the non-abelian case a rest remains that is partly t\ <-> U\ antisymmetric and 
has 1/e 2 and 1/e poles. 

Finally, we note that the gluon self-energy loops in Fig. |^ (i) are zero for massless 
particles (gluons and light quarks) due to Eq. (0). For the massive quark loop one has 
to calculate two-point functions of the type B(k 2 ,m,m) with k\ = and m 7^ 0. Here 
problems occur when naively applying the decomposition procedure, since one would 
divide by k\ for the coefficients B\ and B 2 \. This is the simplest example of the general 
problem that the standard decomposition breaks down whenever projective momenta do 
not exist p2| . Of course, the self-energy integral is simple enough to be calculated directly 
and is given by 



U^(k 2 2 = 0,m^0) = ig 2 s ~5 ab ti- £ J 



(2?r) n (q 2 -m 2 )[(g + k 2 ) 2 - m 2 } 



9 s 2 1 

~~^fiab(k 2 k 2 — fc 2 5 ,W ') — 



,2 



2 / n 

-+JE- ln(4vr) - In 

e \m 2 



.(B7) 



40 



Appendix C: Virtual plus Soft Coefficients 



We list here the polarized coefficients for the virtual plus soft gluon cross section as defined 
in Eq. 



AL QED = [-t x (2t x + u x )/{tu x ) - ui(ti + 2«i)/(titi)] /4 + [-4m 2 s(2t 2 - t lWl + 2u?) 
+ t lUl (5tl + 2t lUl + 5w 2 )] /(4t?u?) + {/?(2m 2 s + 1 2 + 4t lMl + M 2 )/(4t lMl ) 
+ [t 2 u 2 (3t 2 + At x ui + 3u 2 ) + 4m 4 (£ 4 + 2t^ - 8t 2 w 2 + 2t x u\ + m 4 )] /(4t?it?)} £(2) 
- {f3(2m 2 s + t\ + Ahut + uj)/8 + (-24m 4 + 3t 2 + At x Ux + 3u?)/8} /(*iUi) In 2 x 
+ | - [2m 2 tl(6m 4 + 9m 2 tx + 4t 2 ) + ti(2m 2 + ti)(8m 4 + 9m 2 ti + 2t\)ui 



+ 3t(2m 2 + ti) 2 w 2 ] /(4t 2 t 2 «i) - [st\u x {t x + 2ui) + 4m 4 (^ + 3t 2 «i - t x u\ 



+ 2m 2 (t 4 + t\u x - 3t\u\ - 2t x u\ - u\)] /{2f3st\u\) In x In 



-U 



+ | - [ti«i(3ti + 2m) + 4m 6 (3t 2 + U x Ux + 3u 2 ) + m 2 u 2 (15i 2 + 13ti«i + 8u 2 ) 
+ 2mV(12t 2 + 13ti«i + 9u 2 )] /(4txu 2 u 2 ) - [s^u^ti + u x ) - 4m 4 (^ + t\u x 



3ti« 2 - u\) - 2m 2 (t 4 + 2t\u x + 3t 2 w 2 - *iu? - u{)] /(2/?st 2 w 2 ) lnx) In 



-Ml 



77r 



+ { (m 2 s - tiu x ) / {(3txUx) + (3(2m 2 s + t\ + U x ux + w 2 )/(2tiWi) ln(l + x) } In x 

+ {2m 2 t 2 (ui - tx) + t 2 (2t 2 + 2txUx + u\) - 2m 4 (5t 2 + 2t x u x + u 2 )} /(2tf«i) 

x Li 2 (—^J + {2m 2 (t! - u x )u\ + u\{t\ + 2ti«i + 2u\) - 2m A {t\ + 2t x u x + 5u^)} 

/(2tiu?)Ii 2 (^) + p(2m 2 s + tj + Atxu x + u\) / (2t x u x )U 2 (- x) + AB QED 

1 + ln 



x 



1+ i^!(2C(2) + 



— In x 



/3s 

In x + 4Li 2 (x' 



QED 



/2 , 
s — 2m 2 



In x Mi) I — 



_/| ^ +lnf^]+41n(l-x) 



77T 



/5s I \ ;;/ 

2 /j2 , 2\ //r,.2 „.2\ , I /3/<v__2„ , „2 



(CI) 
(C2) 



ALok = m l s{t{ + u{) /{2t{u{) + | - /3(2wTs + + 2t 1 w 1 )/(4£ 1 ui) + ^ - 2tfuJ (2tf 
+ tjuj + 2u^) + m 2 st 1 u x (7t 2 1 - 8t x u x + 7it?) - m A {t\ + 2t?«i - 26t 2 u 2 + 2t x u\ 



+ it 4 ) /(2*?u?) C(2) + (-2m 2 s + t 2 + u 2 x)/{U x ux) 



In 2 



- 1 * + In 2 



-Mi 



- {24m 4 - 3s 2 + 2t x ux - (3{2m 2 s + s 2 + 2ti«i)} /{8t x u x ) In 2 x 



41 



+ <^ (m 2 s + t\)(-m 2 s + t 1 u 1 )/{2U\ui) + \tiUi{t\ + u\) - 2m 2 s{2t\ - t x u x + 2u 2 )] 



/(2tfui)ln 



-Ml 



{st 2 {s - u^m + 2m 2 [st\ + {s 2 + 2t\)u 2 ] - 4m 4 [s 3 + 2t 1 (s' 2 



+ iiwi)]}/(2^ui)lnx|ln ( + |(-m 2 s + tiWi)(m 2 s + w 2 )/(2tiww 2 ) 

- {s(s - ti)tiw 2 - 4m 4 [s 3 + 2 Ul {s 2 + txUx)] + 2m 2 [su\ + t\(s 2 + 2m 2 )] } / (2f3st\ul) 



x In x > In 



x 



-Ui 



■{m 2 s - £iiti)/03ti«i) - /3(2m 2 s + s 2 + 2t 1 w 1 )/(2t 1 « 1 ) 



ln(l + x)} In x + {m 2 t 2 (ti - 3u x ) - tf (ti + 2«i) + m 4 (5t 2 + 2tiiti + w 2 )} /(2tfyi) 
Li 2 f-^J + {m 2 M 2 (-3ti + «i) - wf(2ti + Ml ) + m 4 (t 2 + 2t lMl + 5m 2 )} /(2t lM 3 ) 

x Li 2 - f3{2m 2 s + s 2 + 2t 1 u 1 )/(2t 1 « 1 )Li 2 (-x) + AB QED l - 3((2) 



I In ( 4^1 In (^) + In 2 (± )- 2 In ( h. )\ n x _ \ n 2 x + 2Li 2 (l - 



2Li 2 ( 1 - 



/4 



s — 2m 2 



2<(2) + < ln[-4l+ln 



-Mi 



+ 41n(l — x) — In x > lnx + 4Li 2 (x 



/2 



AL 



OK 



AR 



QED 



"i 



s - 2m 2 , ] , /A 
in x >ln — - 



AL RF = AB QED In 



/4- 



(C3) 
(C4) 
(C5) 



with ABqed given in Eq. ([13]) , /3 = yl — 4m 2 /s, x = (1 — /3)/(l + /3), and t and u as 
denned in (0). 

When integrated, the hard gluon cross section diverges logarithmically ~ ln fc A/m 2 
{k = 1, 2) as the IR cutoff A — > 0. This is by definition canceled by the sum of the L A 
contributions in the soft gluon cross section. If one is interested to show the contributions 
from the hard and the soft plus virtual contributions separately (as, e.g., in Fig. 1 of 
||20||), it is thus advisable to shift the L A terms to the, in this way redefined, hard cross 
section in order to achieve a numerical stable result independent of A. For any numerical 
calculation of physically relevant hadronic cross sections, it is however more practical to 
directly add the complete soft plus virtual piece to the hard cross section. In both cases 
this can be achieved by rewriting the soft plus virtual cross section, expanded in powers 
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of ln fe A/m 2 (k = 0, 1, 2), as follows g|| 

2 



5(s 4 ) y~] a k ln fc — - 9(s 4 - A)A fc a fc | S4=0 (C6) 

k=0 



with certain coefficients at&. (|C6| ) takes properly care of the different distributions <5(s4) 
and 0(s4 — A) multiplying the soft and hard parts, respectively, see Eq. (|34"1). As indicated 
in (|C6j), the a*, have to be always evaluated using the "elastic" 2^2 kinematics, i.e., 
■S4 = 0, even when added to the 2 — > 3 hard cross section. The coefficients Ak are given 
by 



1 ln^f^/m 2 ) 1 ln 2 (s^ ax /m 2 ) 21n(s 4 / 



2^ 



A = r , A 1 = v 4 , A 2 = — v 4 ^ i (C7) 

^max 

as can be easily verified by integrating the l.h.s. of ( |C6|) over (J A 4 ds^), which recovers 
the ln fc A/m 2 terms. 
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